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摘 要 : 为 揭示 呼伦贝尔 沙 地 樟 子 松 人 工 林 土 壤 细菌 相互 关系 ,以 呼伦贝尔 沙 地 不 同 林 龄 樟 子 松 人 工 林 (25 a .34 


a 和 43 a) 为 研究 对 象 ,以 沙 质 草地 为 对 照 ,采用 分 子 生态 网 络 分 析 法 对 不 同 土屋 (0~10 cm 和 10~20 cm ) 土 壤 细 菌 群 


落 进行 比较 分 析 。 结 果 表 明 :(1) 从 25 a 到 43 a, 土 壤 细 菌 网 络 总 边 数 增多 ,平均 路 径 长 度 降低 。 土 壤 深 度 由 0~10 
cm, 人 工 林 土壤 细菌 网 络 总 边 数 减少 ,平均 路 径 长 度 升 高 。 与 沙 质 草地 相 比 , 人工 林 土壤 细菌 网 络 总 
边 数 较 少 。(2) 25 a 人 工 林 关键 菌 种 隶属 于 嗜 酸 菌 日 (Acidimicrobiales) .RB47] 和 MB-A2-108,34 a 人工 林 关键 菌 种 


cm 到 10~20 


隶属 于 Caiellales,43 a 人 工 林 关键 菌 种 隶属 导 


关键 菌 种 隶 
中 介 中 心性 


F Gaiellales , RB41 ,Subgroup. 7 ,Subgroup_6 , #0 DA101_soil_group ,草地 


的 部 分 细菌 有 显著 正 相 关 影 响 (P<0.05) ;转化 酶 和 过 氧化 氢 


属于 匿 杆 菌 门 (Latescibacteria)。(3) ER AR .微生物 碳 含量 和 脲酶 酶 活性 对 土壤 细菌 网 络 中 具有 高 


酶 活性 .土壤 含水 量 以 及 速效 磷 含 量 对 土 


壤 细 菌 网 络 中 具有 高 中 介 中 心性 的 部 分 细菌 有 显著 负 相 关 影响 (P<0.05) ;土壤 有 机 质 对 土壤 细菌 网 络 中 具有 高 中 


介 中 心性 的 部 分 细菌 既 有 显著 正 相 关 影 响 又 有 显著 负 相 关 影响 (P<0.05)。 榜 子 松 人 工 林 从 25 a 到 43 a, ERATA 


网 络 愈加 复杂 和 紧密 ,土壤 深度 由 0~10 cm 到 10~20 cm, 网 络 复杂 性 和 紧密 度 降低 ;与 草地 相 比 ,人 工 林 土 壤 细 菌 网 
络 复杂 性 较 低 。43 a 人 工 林 土壤 细菌 网 络 关键 菌 种 类 型 数量 最 多 。 男 外 ,土壤 细菌 网 络 受 土壤 有 机 质 影 响 最 大 。 
研究 结果 有 助 于 深入 理解 呼伦贝尔 沙 地 榜 子 松 人 工 林 土壤 细菌 群落 ,并 为 呼伦贝尔 沙 地 榜 子 松 信 工 林 的 可 持续 经 


营 提 供 科 技 支 撑 。 


关键 词 : 土壤 微生物 ; 分 子 生 态 网 络 ; 林 龄 ; 土 层 ; 樟 子 松 ; 呼伦贝尔 沙 地 


土壤 细菌 是 土壤 微生物 组 的 重要 类 群 ,也 是 陆 
地 生态 系统 的 重要 组 分 "。 土 壤 细 菌 驱动 陆地 生态 
系统 物质 循环 和 能 量 流动 ,在 生态 系统 稳定 、 调 节 
和 修复 方面 发 挥 着 不 可 符 代 的 重要 作用 ”。 一 般 认 
为 ,气候 .植被 .土壤 等 环境 因素 是 影响 土壤 细菌 群 
落 组 成 和 功能 动态 的 关键 因子 。 此 外 ,土壤 细菌 
间 的 相互 关系 与 细菌 群落 结构 也 蜜 不可分。 研究 
表明 ,土壤 细 戎 群落 结构 庞大 ,相互 关系 复杂 ,土壤 
细菌 之 间 既 有 积极 的 协同 作用 ,又 有 消极 的 持 抗 作 
用 ,这 些 复杂 的 相互 关系 在 细菌 群 落 的 塑造 上 发 挥 
着 重要 作用 ,同时 对 维持 生态 系统 功能 和 服务 的 影 
响 也 十 分 重要 1。 但 是 ,由 于 传统 土壤 细菌 相互 关 
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系 识别 和 判定 方法 的 局 限 性 ,使 我 们 对 于 土壤 细菌 
相互 关系 依然 缺乏 全 面 深入 的 了 解 。 

分 子 生态 网 络 分 析 是 一 种 新 兴 的 探索 性 数据 
分 析 方 法 ,该 方法 基于 随机 和 矩阵 理论 ,结合 基因 测 
序 技术 ,将 微生物 间 的 交互 作用 以 及 微生物 系统 稳 
定性 可 视 化 。 分 子 生态 网 络 分 析 技 术 成 为 解析 微 
生物 间 相 互 关 系 及 其 影响 因素 的 重要 手段 " ,该 方 
法 在 净值 控制 下 可 揭示 微生物 组 中 复杂 的 戎 群 关 
系 以 及 鉴别 关键 菌 种 。 研 究 发 现 , 中 国 北方 森林 
较 南方 森林 土壤 微生物 群落 间 关 系 更 复杂 ,南方 
红壤 区 马尾 松林 土壤 微生物 群落 以 负 相 关 关 系 为 
主 ,中国 北方 沙漠 土壤 细菌 和 真 核 生物 相互 作 
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用 , 随 着 年 平均 降水 量 的 增加 而 增强 ”。 细 菌 网 络 内 
特定 的 相互 作用 可 能 增加 或 减少 生物 多 样 性 ”” ,也 
正 是 因为 具有 这 些 相 互 作用 ,细菌 网 络 才 可 以 更 好 
的 调度 各 种 土壤 细菌 ; 当 群 落 受 到 扰动 时 ,细菌 网 
络 可 以 招募 具有 不 同 细菌 组 合 的 替代 途径 来 实现 
特定 功能 ,具有 缓冲 作用 ”。 因 此 ,分 子 生态 网 络 分 
析 是 进一步 解析 土壤 细菌 之 间 相 互 作用 及 其 影响 
因素 的 重要 手段 ,也 是 揭示 土壤 细菌 网 络 结构 动态 
与 生态 过 程 和 作用 协同 变化 的 重要 途径 。 

ta TW (Pinus sylvestris var. mongolica) 是 我 国 北 
方 风沙 区 防风 固沙 的 重要 树种 ,土壤 细菌 间 相 互 关 
系 对 樟 子 松林 木 生长 和 健康 状态 具有 重要 影响 “” , 土 
壤 细 菌 网 络 的 “缓冲 "作用 ,在 榜 子 松林 抵抗 环境 变 
化 方面 发 挥 着 重要 作用 "。 鉴 于 此 ,以 呼伦贝尔 沙 
地 25 年 生 、34 年 生 和 43 年 生 榜 子 松 入 工 林 为 研究 
对 象 (简称 25 a、34 a 和 43 a) ,以 临近 地 段 沙 质 草地 
为 对 照 , 利 用 分 子 生态 网 络 分 析 技 术 ,揭示 樟 子 松 
人 工 林 土壤 细菌 网 络 结构 特征 及 其 影响 因素 ,以 期 
深入 理解 呼伦贝尔 沙 地 榜 子 松 入 工 林 土 壤 细 菌 相 
互 关 系 ,并 为 区 域 樟 子 松 人 工 林 可 持续 经 营 提供 理 
论 依据 和 科技 支撑 。 


1 材料 与 方法 


1.1 研究 区 概况 

人 研究 区 位 于 内 蒙古 海拉尔 樟 子 松 国 家 森林 公 
be] (119*31' -119?43' E, 49*05' -49*13' N) ,海拔 606 
m, 属 北 温带 大 陆 性 半 干 旱 气 候 。 年 均 气 温 -0.3 C, 
年 均 降水 量 359.4 mm ,年 均 蒸发 量 1148.5 mm, 年 均 
日 照 时 数 2559 h ,无 霜 期 121 d 左 右 。 人 研究 区 主要 十 
壤 类 型 为 暗 栗 钙 土地 带 风 沙土 ,主要 乔木 植被 有 樟 
FHS LEER BK (Quercus mongolica ) 、 杨 树 (Populus to- 
mentosa ) , AE (Betula platyphylla) 等 ,主要 草本 植物 
A E We 33€ (Potentilla chinensis) , £F $ (Stipa capilla- 


ta) . > (Leymus chinensis ) 、 狗 尾 草 (Setaria viridis ) 
Fibi (Sanguisorba officinalis ) 5 5 
12 样品 采集 与 处 理 

2018 年 7 月 ,在 内 蒙古 海拉尔 樟 子 松 国 家 森林 
公园 选取 无 经 营 管理 措施 、 林 下 植被 相似 的 25 a、 
34 a 和 43 a 樟 子 松 人 工 林 为 研究 对 象 ,同时 选取 附 
近 沙 质 草地 作为 对 照 。 在 不 同 林 龄 的 樟 子 松 人 工 
林 及 沙 质 草地 内 分 别 随 机 布设 50 mx50 m 的 样 地 
( 表 1)。 样 地 内 进行 每 木 检 尺 后 每 块 樟 子 松 人 工 林 
样 地 随机 选取 3 株 标准 木 ,在 树冠 投影 处 采集 土壤 
样品 ,采集 深度 为 0~10 cm 和 10~20 em; 沙 质 草地 样 
地 随机 选取 3 个 取样 点 采集 土壤 样品 ,样品 采集 次 
BEJ 0~10 cm 和 10~20 em。 每 个 取样 点 的 同一 土 层 
随机 采集 3 个 样品 后 充分 混合 为 1 个 样品 ,每 份 土 
壤 样 品 约 为 500 g, 共 24 份 土壤 样品 。 采 样 期 间 为 
E 1E ERRE mIa e , ECT IURE H 7590 URB JC TRI 
消毒 。 同 时 ,分 别 用 环 刀 和 铝 盒 采集 土壤 样品 ,用 于 
土壤 含水 量 的 测定 。 所 有 土壤 样品 置 于 4 %C 便 携 式 
保温 箱 保存 , 运 回 实验 室 后 转 和 -80 冰箱 中 保存 。 
1.3 土壤 理化 性 质 测定 

采用 烘 干 法 测定 土壤 含水 量 (Soil Water Con- 
tent,SWC)。 其 余 土壤 样品 风干 后 过 2 mm ffi, HF 
土壤 样品 化 学 性 质 分 析 : 采 用 靛 酚 蓝 比 色 法 测定 全 
ZA (Total Nitrogen, TN) M A A (Ammonia Nitrogen, 
AN) KJ CHR 2 FF i PEA DU E ERA HLE (Soil 
Organic Matter, SOM) , FH 4H ERDE EG IS I E UR 
BE ( Available Phosphorus, AP), 采用 高 锰 酸 钾 滴 定 法 
测定 过 氧化 氧 酶 含量 (Catalase Activity, CA) ,采用 
茶 酚 钠 - 次 握 酸 钠 比 色 法 测定 脲酶 活性 (Soil Ure- 
ase, SU) ,采用 NasS;0; 滴 定 法 测定 转化 酶 活性 (In- 
vertase Activity, LA) ,采用 氯仿 束 菊 法 测定 微生物 碳 
(Microbial Carbon, MC) ^^, 
1.4 土壤 细菌 的 分 子 鉴 定 

使 用 PowerSoil DNA Isolation Kit 试剂 盒 (MoBio 


R1 呼伦贝尔 沙 地 模子 松 信 工 林 样 地 概况 


Tab.1 Basic information of P sylvestris plantations sample plots in the Hulunbuir Desert 


样 地 编号 林 龄 /a 平均 树 高 /m 平均 胸径 /cm 林 分 密度 /( 株 hm?) HBAJ RE 土壤 类 型 
HLP25 25 12.16+1.21 15.24+2.79 1650 0.76 风沙 土 
HLP34 34 12.71+0.97 19.25+3.10 1650 0.84 风沙 土 
HLP43 43 13.52+1.24 23.59+3.28 1650 0.73 风沙 土 

HLG - - - - 0.75 风沙 土 


注 :HL 表 示 呼 伦 贝尔 ;G 表示 沙 质 草地 ;P 表 示人 工 林 。 下 同 。 
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Laboratories, Carlsbad, CA ) 对 土壤 细菌 样品 进行 
DNA 抽 提 ,完成 基因 组 DNA 抽 提 后 ,利用 1% 琼 脂 糖 
凝 胶 电泳 检测 抽 提 的 基因 组 DNA。 随 后 进行 PCR 
扩 增 , 扩 增 区 域 为 16S rRNA 基因 V3-V4 区 ; 扩 增 引 
物 序 列 为 338F (5'CCTACGGGAGGCAGCAG- 3’) , 
806R (5'-ATTACCGCGGCTGCTGG-3) ; 2 3E A EA 
95 CHAE TE 300 s, 95 CARTE 45 s,55 "OB K 50s, 
72 CHEM 45 s,28 个 循环 ,土壤 细菌 样品 分 两 层 取 
样 ,每 层 土壤 细菌 样品 做 3 次 重复 ,将 同一 样品 的 
DNA 混合 ,作为 每 份 土 壤 细 菌 样品 的 总 DNA。 将 土 
壤 细 菌 样品 的 PCR 产物 混合 后 使 用 2% 琼 脂 糖 凝 觅 
电泳 检测 ,使 用 AxyPrepDNA 凝 胶 回收 试剂 盒 (AX- 
YGEN Zt 8] ) Je dit PCR 7^9 , FEF Tris-HCl Z vp 
T UE, Tee fi FH 296 BA Fl, v e I 

Illumina Miseq X 库 构 建 之 后 进行 lumina 
Miseq 上 机 测序 ,测序 完成 后 使 用 Trimmomatic、 
FLASH , Pear, usearch 软件 对 数据 进行 处 理 , 基 于 
Silva (Release128/132 http://www.arb- silva.de) , RDP 
(V16 http://rdp.cme.msu.edu/) 和 Greengene (Re- 
leasel3.5 http://greengenes.secondgenome.com/ ) 数据 
JE ,采用 RDP Classifier $5 Y% BY blast, uclust consen- 
sus taxonomy assigner 等 方法 对 OTU 代表 序列 进行 
比 对 分 析 ,并 对 OTU 聚 类 统计 的 相似 水 平 控制 在 
97% ,得 到 分 类 群 的 生物 信息 “。 测 序数 据 NCBI 登 
录 号 为 PRJNA875585。 
1.5 数据 分 析 

筛选 相对 丰 度 大 于 5.00% 的 菌 属 ,保留 其 中 在 
每 个 样本 中 都 出 现 的 OTU ,使 用 SPSS 24.0 计 算 样品 
OTU FY) Pearson 相关 性 与 显著 性 值 (r>0.6 和 P<0.05)， 
之 后 使 用 RStudio 计算 Pearson FAK ABE. fE 
HH RStudio*igraph ”" 包 将 相关 系数 和 矩阵 转化 为 网 络 文 
件 并 计算 其 总 边 数 总 点 数 .平均 聚 类 系数 .平均 路 
BRE .模块 化 .平均 连通 度 .中 介 中 心性 等 拓扑 特 
征 ,同时 计算 模块 内 连通 度 (2) 和 模块 间 连 通 度 
(Pi) 并 进行 网 络 节点 的 模块 角色 判别 ,包括 模块 
中 心 点 (module hubs) 、 网 络 中心 点 (network hubs ) , 
连接 节点 (connectors) 和 外 围 节 点 (peripherals)。 通 
常 ,Zi2.5 或 Pi>0.62 的 节点 被 认为 是 关键 物种 ( 表 
2) ,其 在 模块 间 和 各 自 的 模块 内 起 重要 作用 。 使 用 
Gephi 0.9.2 将 网 络 文件 可 视 化 并 导出 土壤 细菌 共 现 
网 络 图 ,网 络 图 中 的 每 一 个 节点 (node) 代 表 一 个 
OTU ,连接 市 点 之 间 的 每 条 边 (edge) 代 表 各 0TU 之 
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表 2 网 络 节 点 的 判别 


Tab.2 Discrimination of network nodes 


网 络 节 点 模块 内 连通 度 (Zi) 模块 间 连 通 度 (Pi) 
模块 中 心 点 >2.5 <0.62 

连接 节点 «2.5 20.62 
网 络 中 心 点 >2.5 >0.62 

外 围 节 点 <2.5 <0.62 


间 的 相关 性 ,基于 各 0TU 之 间 的 相关 性 ,形成 相互 
联系 的 模块 ,并 计算 模块 化 (modularity ) 数 值 (模块 
化 大 于 0.44 时 网 络 图 有 意义 )。 最 后 利用 RStudio 
“igraph” 包 生成 1000 个 节点 和 边 数 与 共 现 网 络 相同 
的 随机 网 络 ,计算 拓扑 参数 ,比较 其 与 细菌 共 现 网 
络 的 平均 聚 类 系数 ,判别 土壤 细菌 网 络 凝 聚 性 是 否 
有 意义 。 

为 了 分 析 细 菌 生态 网 络 与 土壤 因子 的 关系 , 选 
择 樟 子 松 人 工 林 土壤 细菌 网 络 最 具 代 表 性 的 OTU 
(中 介 中 心性 >9000 的 OTU ) 与 土壤 因子 建立 网 络 关 
系 。 以 中 介 中 心性 值 与 土壤 因子 进行 Person 相关 性 
分 析 (r>0.6 和 P<0.05) 并 使 用 Gephi 0.9.2 将 土壤 细 
菌 与 土壤 因子 互 作 图 可 视 化 ,最 后 导出 土壤 细菌 与 
土壤 因子 相互 作用 图 。 


2 结果 与 分 析 


2.1 樟 子 松 人 工 林 土壤 细菌 网 络 高 连通 性 的 拓扑 
性 质 


与 


呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 细菌 群落 经 
验 网 络 模块 化 .平均 路 径 长 度 和 平均 聚 类 系数 明显 
高 于 随机 网 络 ,符合 网 络 特征 ,可 用 于 土壤 细菌 互 
作 的 研究 ( 表 3)。 

樟 子 松 信 工 林 土 壤 细 菌 群 落网 络 拓扑 特征 存 
在 明显 林 龄 差异 ( 表 3, 图 1)。 从 25 a 到 43 a, 土 壤 
细菌 网 络 总 边 数 .平均 连通 度 增 加 ,总结 点 数 和 平 
均 路 径 长 度 降 低 ; 相 比 之 下 ,草地 土壤 细菌 网 络 总 
结 点 数 和 总 边 数 多 于 人 工 林 , 但 平均 路 径 长 度 比 人 
工 林 要 短 。 总 边 数 越 多 ,网 络 愈加 复杂 ;平均 路 
径 长 度 越 低 , 网 络 愈加 紧密 。 从 25 a 到 43 a, 樟 子 松 
人 工 林 土壤 细菌 网 络 愈加 复杂 和 紧密 ;与 人 工 林 相 
比 ,草地 土壤 细菌 网 络 更 加 复杂 和 紧密 。 

樟 子 松 信 工 林 土 壤 细 菌 群 落网 络 拓扑 特征 存 
在 明显 土 层 差异 ( 表 3, 图 1) ,表层 土壤 细菌 网 络 较 
下 层 总 边 数 平均 连通 度 、 正 相关 关系 占 比 更 高 , 平 
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表 3 呼伦贝尔 沙 地 樟 子 松 信 工 林 土壤 细菌 网 络 拓 扑 特征 


Tab.3 Topological characteristics of soil bacterial network of P sylvestris plantations in the Hulunbuir Desert 


拓扑 特征 HLP25A HLP25B HLP34A HLP34B HLP43A HLP43B HLGA HLGB 
经 验 网 络 节点 306 304 307 292 300 268 330 325 
边 2395 2262 2483 2313 3072 2011 2689 2831 
平均 连通 度 15.65 14.88 16.18 15.84 20.48 15.01 16.30 17.42 
模块 化 0.75 0.77 0.70 0.72 0.62 0.62 0.71 0.71 
平均 路 径 长 度 11.89 15.26 11.50 11.67 8.57 13.81 12.14 11.52 
平均 聚 类 系数 0.75 0.74 0.75 0.75 0.76 0.74 0.74 0.75 
正 相 关 /% 54.53 51.50 52.56 49.72 52.96 50.22 48.83 53.27 
负 相 关 /% 45.47 48.50 47.44 50.28 47.04 49.78 51.17 46.73 
随机 网 络 模块 化 0.21 0.21 0.20 0.20 0.18 0.21 0.21 0.20 
平均 路 径 长 度 2.37 2.41 2.35 2.34 2.16 2.35 2.37 2.32 
平均 聚 类 系数 0.05 0.05 0.06 0.05 0.07 0.05 0.05 0.05 


注 :A 表 示 0~10 cm 土 层 ;B 表 示 10~20 cm 土 层 。 下 同 。 


HLP25B 
注 :A 表 示 0~10 cm 土 层 ;B 表 示 10~20 cm 
性 ,绿色 为 正 相 关 , 红 色 为 负 相 关 。 


HLP34B 


FT 


HLP43B 


HLGB 


E. PAPAS ET es BUG AR BEES [nd ; 节点 大 小 表示 中 介 中 心性 大 小 ; 边 的 颜色 表示 正 负 相关 


图 1 呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 细菌 共 现 网 络 


Fig. 1 Co-occurrence network of soil bacterial communities of P. sylvestris plantations in the Hulunbuir Desert 


均 路 径 长 度 更 低 。 与 之 相反 ,草地 表层 土壤 细菌 网 
络 较 下 层 总 边 数 平均 连通 度 、 正 相关 关系 占 比 更 
低 , 平 均 路 径 长 度 更 高 。 土 壤 深 度 由 0~10 cm 到 10~ 
20 cm, 人工 林 土 壤 细 菌 网 络 复杂 性 和 紧密 度 降 低 ; 
草地 土壤 细菌 网 络 复杂 性 和 紧密 度 增 加 。 
2.2 樟 子 松 土壤 细菌 网 络 的 模块 特性 和 关键 菌 种 
樟 子 松 人 工 林 和 草地 土壤 细菌 共 现 网 络 外 转 
节点 占 总 节点 的 99% 以 上 ,未 检测 到 模块 中 心 点 和 


a 


网 络 中 心 点 ,共识 别 到 11 个 连接 节点 。 其 中 ,人 工 
林 连 接 节点 10 个, 分别 隶属 于 MB-A2-108、Sub- 
group_6、 嗜 酸 菌 日 (Acidimicrobiales) 、Gaiellales、 
Subgroup_7、DA101_soil_group 和 RB47; 草 地 连接 节 
点 1 个 ,隶属 于 匿 杆菌 门 ( 表 4, 图 2)。 
2.3 樟 子 松 人 工 林 土壤 细菌 网 络 的 影响 因素 
呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 理化 性 质 和 
酶 活性 波动 性 较 强 ( 表 5)。 从 25 a 到 43 a, 表 层 土壤 
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表 4 呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 细菌 网 络 连 接 节 点 


Tab.4 Connectors of soil bacterial network of P. sylvestris plantations in the Hulunbuir Desert 


关键 菌 种 OTU 门 2 H 科 属 
HLP25A OTU144 Actinobacteria Acidimicrobiia Acidimicrobiales = - 
HEPIE OTU182 Acidobacteria Blastocatellia Blastocatellales Blastocatellaceae Subgroup. 4 RB4I 

OTU291 Actinobacteria MB-A2-108 = = B 
HLP34A OTU149 Actinobacteria Thermoleophilia Gaiellales - - 
HLP34B - - - - - - 
OTU66 Acidobacteria Blastocatellia Blastocatellales Blastocatellaceae_Subgroup_4 RB4I 
OTU677 Acidobacteria Holophagae Subgroup. 7 i = 
HLP43A 
OTU319 Acidobacteria Subgroup_6 = = - 
OTUS556 Verrucomicrobia Spartobacteria Chthoniobacterales DAIOIL soil group - 
HLPA3B OTUI52 Acidobacteria Subgroup. 6 = = = 
OTU3844 Acidobacteria Thermoleophilia Gaiellales - - 
HLGA - - - - - - 
HLGB OTU718 Latescibacteria - ka - = 
SERA A HURE et w T EE CP«0.05) ,微生物 碳 含 T 
量 先 升 高 后 降低 (P<0.05) ,过 氧化 氢 酶 和 脲酶 活性 ”3 讨 E 


先 降低 后 升 高 (P<0.05 )。 沙 质 草 地 土壤 全 和 毛 含 量 显 
著 高 于 人 工 林 (P<0.05) ,脲酶 活性 显著 低 于 人 工 林 
(P<0.05)。 从 25 a 到 43 a, 下 层 土壤 全 所 和 微生物 
碳 含量 先 上 升 后 下 降 (P<0.05) ,氨氮 含量 显著 下 降 
(CP<0.05)。 沙 质 草 地 土壤 转化 酶 活性 显著 高 于 人 工 
林 (P<0.05)( 表 5)。 

呼伦贝尔 沙 地 土壤 细菌 网 络 中 具有 高 中 介 中 
心性 的 部 分 菌 种 与 土壤 含水 量 ` 有 机 质 .速效 磷 、 氮 
AAA 微生物 碳 、 氧 化 氢 酶 .转化 酶 和 脲酶 显著 
相关 (P<0.05) ,其 中 土壤 有 机 质 影响 最 大 (图 3)。 
土壤 有 机 质 与 微 枝 形 杆菌 属 (MWicrovirga) 的 OTU235 
和 Burkholderia- Paraburkholderia 的 OTU451 =Œ fj E 
正 相 关 (P<0.05) , 5j Gaiella 的 OTU29 Fil Acidobacte- 
riaceae, Subgroup. 1 的 OTU544 Œ fij # f 4H X (P< 
0.05) ;土壤 含水 量 和 转化 酶 与 芽 单 胞 菌 科 (Gemma- 
timonadaceae ) I OTU145 呈 显 著 负 相关 (P<0.05 ) ; 4& 
氮 与 RB41 的 0TU256 呈 显著 正 相 关 (P<0.05 ) ; BA 
与 Subgroup_6 的 OTU194 和 DA101_soil_group 的 
OTU2222 呈 显 著 正 相关 (P<0.05) ; 速效 磷 与 Acido- 
bacteriaceae, Subgroup. 1 的 OTU155 和 MB-A2- 108 
的 OTU432 E tj E ff HOS (P<0.05) ;过 氧化 氧 酶 与 
Subgroup_6 的 OTU12 星 显著 负 相 关 (P<0.05) ;脲酶 
与 Burkholderia- Paraburkholderia Ef] OTU451 Œ ii 3E 
正 相 关 (P<0.05) ;微生物 碳 与 Roseiflexus 的 OTU773 
呈 显 著 正 相关 (P<0.05)( 图 3)。 


3.1 樟 子 松 人 工 林 土 壤 细 菌 网 络 结构 的 特征 

呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 细菌 网 络 拓 
扑 性 质 随 林 龄 呈现 动态 变化 。 从 25 a 到 43 a, 土 壤 
细菌 虽然 种 类 减少 ,但 其 网 络 愈加 复杂 和 紧密 。 研 
究 发 现 ,平均 路 径 长 度 越 短 ,土壤 细菌 网 络 物种 间 
物质 循环 .能量 流动 和 信息 传递 的 效率 越 高 中 ,说 
明 当 环境 受到 较 大 扰动 时 ,人 工 林 土壤 细菌 的 响应 
速度 从 25 a 到 43 a 变 快 , 抗 干 扰 能 力 变 弱 , 细 菌 群落 
稳定 性 变 弱 ""。 这 可 能 与 模子 松 人 工 林 的 衰老 有 
ox", 土壤 细 菌 网 络 愈加 复杂 ,细菌 间 联 系 频繁 可 
能 就 是 为 了 抵抗 人 工 林 衰老 带 来 的 负面 影响 。 

与 樟 子 松 人 工 林 相 比 ,草地 土壤 细菌 种 类 更 为 
丰富 、 相 互联 系 增 多 ,使 得 其 网 络 规模 较 大 且 更 加 
复杂 ,这 与 前 人 研究 结果 一 致 加。 这 是 因为 樟 子 松 
根系 可 以 通过 分 泌 代 谢 物质 招募 对 其 生长 有 利 的 
土壤 细菌 ,排斥 对 其 生长 有 害 的 土壤 细菌 ,宿主 效 
应 种 对 土壤 细菌 进行 筛选 和 过 滤 ,使 得 人 工 林 下 十 
壤 细 菌 多 样 性 和 复杂 性 降低 ”。 男 外 ,草地 土壤 细 
菌 网 络 比 人 工 林 更 加 紧密 ,意味 着 草地 细菌 群落 稳 
定性 比 人 工 林 差 ,这 与 前 人 研究 结果 一 致 ”> 。 

呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 细菌 网 络 在 
不 同 土壤 深度 中 的 拓扑 特征 不 同 。 随 着 土壤 深度 
的 增加 人 工 林 土 壤 细 菌 相互 联系 减少 .平均 路 径 长 
度 增 加 ,网 络 复杂 性 和 紧密 性 降低 ,这 与 前 人 研究 
结果 一 致 站 。 这 是 因为 研究 样 地 0~10 cm 土 层 土壤 
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图 2 呼伦贝尔 沙 地 模子 松 信 工 林 土壤 细菌 拓扑 角色 


Fig. 2 Soil bacterial topological roles of P sylvestris plantations in the Hulunbuir Desert 


细菌 群落 多 样 性 指数 要 高 于 10~20 cm 土 层 ” ,土壤 
细菌 群落 具有 明显 表 聚 性 ”。 因 此 , 相 较 于 0~10 
em 土 层 ,10~20 em 土 层 土壤 细菌 数量 少 .迁移 速率 
AI ,它们 接触 产生 联系 的 机 会 变 少 ,所 形成 的 网 络 
较 玻 松 。 同 时 ,0~10 cm 土 层 营养 物质 丰富 .土壤 孔 
际 较 大 、 含 氧 量 较 高 ,这 些 环境 因素 适合 土壤 细菌 
的 生长 与 交流 “。 男 外 , 共 现 网 络 中 正 负 相关 性 分 
别 可 表示 菌 群 之 间 的 协同 作用 和 持 抗 作用 S E 


二 细菌 间 协 同 作 用 有 助 于 土壤 细菌 生存 繁衍 和 合 
作 分 解 枯 落 物 … ,虽然 土壤 深度 由 0~10 cm 到 10- 
20 cm, 榜 子 松 人 人工 林 土 壤 细 菌 间 协同 作用 逐渐 减 
弱 , 持 抗 作用 增强 ,但 总 体 上 还 是 协同 作用 占 主导 ， 
这 有 助 于 土壤 细菌 群落 更 有 效 地 分 解 枯 落 物 并 为 
榜 子 松 的 生长 提供 养分 。 
3.2. 樟 子 松 人 工 林 土壤 细菌 网 络 关键 菌 种 生态 功能 
不 同 林 龄 土壤 细菌 网 络 关 键 菌 种 不 同 ,存在 物 
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表 5 呼伦贝尔 沙 地 樟 子 松 人 工 林 土壤 理化 性 质 与 酶 活性 


Tab.5 Soil physical and chemical properties and soil enzyme activities of P. sylvestris plantations in the Hulunbuir Desert 


土壤 因子 HLP25A HLP25B HLP34A HLP34B HLP43A HLP43B HLGA HLGB 
土壤 含水 量 SWC/% 10.4423.66a 5.77+1.39A 9.00+1.27a 6.66+0.78A 6.88+0.47b — 5.8122.08A 10.18+0.89a 7.69+1.35A 
全 气 TN/(gkg') 0.82+0.01b C0.56+0.01C 0.77+0.01le 0.63+0.00A — 0.5140.03d 0.49+0.01D 0.93+0.02a 0.59+0.01B 
土壤 有 机 质 SOM 2.31+0.06a 1.90+0.05A 2.11+0.05b 1.82+0.02AB 1.71+0.12c¢ 1.32+0.06C 2.01+0.08b 1.74+0.06B 
/(g+ke"') 
dA AN/(mg*kg') 45.3743.08ab ”84.73+1.59A — 43.3042.19b 44.57+2.99B 50.0643.36a 30.63+5.17C  49.8243.67a 44.12+0.57B 
速效 磷 AP/(g:kg') 3.87+0.3la ——3.35€0.51A — 6.11&1.44a — 4.9340.30A 5.40+1.96a — 5.0740.84A 4.13+0.65a 3.28+1.86A 


过 氧化 氨 酶 CA/(U. 938.76+31.28b545.09+528.63AB 132.00+ 172.64+7.02B 997.9544.60a 921.19+29.80A 992.65+22.58a884.29+25.03A 

g') 38.1 1e 

THENEIA/CU-g") 70.92+6.00c 41.24+2.36B 69.63+1.1lc 41.64+7.94B 97.93+3.00a 48.94+6.24B 81.17+4.79b — 74.6148.92A 

脲酶 SU/(U'g') 635.42+35.36a 213.98+68.25A 396.04+ 251.52+20.26A511.12+26.12b 223.20+ 155.32+28.05d258.03+26.61A 
43.47c 162.05A 

微生物 碳 MC 21.58+12.89c 18.52+14.02B 73.09+1.78a 44.89+2.15A 14.32+1.74c 11.60+3.30B 59.32+4.63b 51.75+1.16A 

/(mg:kg ') 


注 :不 同 小 写 .大 写字 母 分 别 表示 0~10 em £l 10-20 em 土 层 不 同样 地 土壤 细菌 相对 丰 度 差异 显著 (P<0.05)。 


OTUI155 


e OTU194 
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注 : 边 的 颜色 表示 正 负 相关 性 ,绿色 为 正 相 关 , 红 色 为 负 相关 。 
图 3 呼伦贝尔 沙 地 樟 子 松 人 工 林 土 壤 细菌 与 土壤 因子 
网 络 关系 


Fig. 3 Network of soil bacterial and soil factors of 


P. sylvestris plantations in the Hulunbuir Desert 


种 角色 转化 现象 。 不 同 林 龄 阶段 的 土壤 细菌 网 络 
关键 菌 种 存在 差异 原因 可 能 是 土壤 生境 改变 以 及 
宿主 根系 分 泌 物 招募 策略 变化 导致 群落 网 络 构建 
的 动态 变化 。25 a 和 34 a 人 工 林 关 键 菌 种 类 型 较 少 
且 多 隶属 于 放 线 菌 门 SEIS T RB41.MB-A2-108, 
Gaiellales 和 嗜 酸 菌 目 的 部 分 细菌 在 樟 子 松 与 土壤 
碳 氮 磷 交 流 中 发 挥 重要 作用 。 其 中 隶属 于 嗜 酸 
菌 目 的 细菌 能 够 将 植物 天 然 产生 的 且 分 解 为 羧 酸 
盐 和 锭 最 后 以 碳 源 和 氮 源 的 形式 被 微生物 和 植物 
AFA), xx nIBEREDSI 25 a 和 34 a 人 工 林 分 别处 于 


中 林 龄 和 近 熟 林 阶 段 ,人 工 林 与 土壤 细菌 间 养 分 交 
流 密切 ,有 利于 促进 植物 生长 ””。43 a 人 工 林 关 键 
AMKA RE ,主要 隶属 于 酸 杆菌 门 , 除 隶 属于 Gai- 
ellales 和 RB41 的 OTU 外 ,43 a 人 工 林 关 键 菌 种 还 检 
测 出 隶属 于 Subgroup 6, Subgroup. 7、 和 DA101_ 
soil. group 的 OTU。43 a 人工 林 处 于 成 熟 林 阶 段 , 养 
分 归还 成 为 重要 生态 过 程 ,隶属 于 酸 杆菌 门 的 关键 
菌 种 在 植物 纤维 素 的 分 解 上 扮演 着 重要 角色 。 
男 外 ,有 人 研究 表明 DA101_soil_group 是 干旱 草地 的 
优势 菌 种 并 且 干 旱 胁 迫 会 导致 Subgroup_7 P8 ES] 
Ee, 

不 同 植被 土壤 细菌 网 络 关 键 菌 种 类 型 存在 差 
异 。 人 工 林 土壤 细菌 网 络 关 键 菌 种 在 类 型 上 优 于 
草地 。 关 键 阔 种 类 型 越 多 ,网 络 稳定 性 越 高 ”。 从 
关键 菌 种 角度 出 发 ,人 工 林 网 络 稳定 性 高 于 草地 ， 
这 与 本 文 研 究 结 论 一 致 。 草 地 土壤 细菌 网 络 关键 
菌 种 隶属 于 匿 杆 菌 门 , 匿 杆菌 门 细 菌 普遍 存在 腐生 
的 生活 方式 , 且 偏 好 大 氧 和 富 营 养 化 的 栖息 地 。 
这 说 明 草 地 土壤 含 氧 量 和 和 孔 阶 度 可 能 要 低 于 人 
工 林 。 
3.3 土壤 因子 对 土壤 细菌 网 络 的 调控 作用 

土壤 理化 性 质 是 呼伦贝尔 沙 地 土壤 细菌 群落 
的 驱动 因素 之 一 。 其 主要 受 土壤 有 机 质 调控 。 研 
究 表明 ,植被 很 大 程度 上 影响 土壤 性 质 ,尤其 是 土 
壤 碳 含量 "。 有 机 质 作为 土壤 养分 库 的 重要 组 成 
部 分 ,与 不 同 土壤 细菌 时 现 不 同 的 相关 关系 ,是 
影响 土壤 细菌 群落 结构 的 主要 因素 ” ,这 与 本 研究 
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结果 一 致 。 土 壤 有 机 质 的 增加 和 减少 都 可 刺激 土 
壤 细 菌 的 生命 活动 ,土壤 细菌 通过 一 系列 的 生理 反 
应 影响 着 植物 对 所 和 磷 的 吸收 和 利用 ,土壤 有 机 质 
的 供应 对 和 森林 的 生长 发 育 至 关 重 要 。 随 林 龄 增 
长 , 樟 子 松 吸收 和 固定 碳 的 能 力 、 枯 落 物 养分 归还 
砚 的 能 力 有 着 显著 差异 ,土壤 细菌 群落 主要 从 土 
二 有 机 质 分 解 中 获得 碳 源 和 能 量 ,其 活性 主要 受 土 
壤 碳 含量 和 分 解 过 程 中 微生物 碳 含量 的 影响 ,土壤 
细菌 群落 在 樟 子 松林 生长 的 每 个 阶段 都 在 通过 调 
整 群 落 内 部 细菌 活性 来 适应 变化 的 环境 ”。 在 此 
过 程 中 ,部 分 土壤 细菌 会 断 开 或 是 重新 建立 与 其 他 
细菌 的 联系 ,这 是 一 种 更 好 的 生存 策略 ”"。 速 效 磷 
含量 的 变化 显著 影响 在 土壤 细菌 网 络 中 发 挥 重要 
功能 的 某 些 菌 种 的 同时 ,氨氮 和 全 气 对 土壤 细菌 
网 络 也 具有 正 向 调控 作用 ,这 是 因为 氮 、 克 元 素 由 
是 陆地 生态 系统 最 主要 的 限制 性 营养 元 素 又 是 能 
够 影响 土壤 细菌 群落 代谢 活性 的 重要 元 素 "“。 土 
SEE Dre ABETE Ae AY VE E ht FAIL AS RSC) AE 
PRAY BE 22H UA , ERWA ER He B AE f n] 
使 土壤 细菌 互 作 种 类 发 生变 化 (如 持 抗 和 协同 作用 
的 转变 ) ,从 而 改变 土壤 细菌 网 络 复杂 性 A 
外 , 氮 须 和 有 机 质 含 量变 化 对 微生物 碳 含量 有 重要 
影响 ” ,微生物 碳 可 反映 土壤 生物 活性 ,也 是 土壤 
细菌 网 络 的 支撑 基础 。 全 毛 、 氨 所 和 速效 磷 分 别 与 
隶属 于 RB41、Subgroup_6 和 MB-A2-108 的 细菌 有 
显著 相关 性 (P<0.05) ,而 这 些 细菌 的 分 类 水 平 与 关 
刍 菌 种 的 分 类 水 平 有 重 闭 ,他 们 可 能 有 着 相同 或 类 
似 的 生态 作用 。 当 外 界 环 境 发 生 扰 动 时 ,土壤 细菌 
网 络 可 以 招募 这 些 细 菌 代替 关键 菌 种 发 挥 特定 生 
态 功能 ,以 此 来 达到 缓冲 作用 。 

酶 也 是 呼伦贝尔 沙 地 土壤 细菌 群落 的 驱动 
素 之 一 。 土 壤 细 菌 群落 调控 因子 的 多 样 性 预示 着 
呼伦贝尔 沙 地 生态 平衡 不 易 被 打破 ”。 随 着 沙 地 土 
培养 分 的 变化 ,细菌 代谢 途径 也 变 得 多 样 。 细 戎 可 
分 泌 多 种 酶 来 降解 这 些 复 杂 多 样 的 营养 物质 。 本 
研究 指出 脲酶 .过 氧化 氢 酶 和 转化 酶 是 呼伦贝尔 沙 
地 物质 循环 中 不 可 或 缺 的 重要 酶 种 ,扮演 着 重要 角 
色 。 呼 伦 贝尔 沙 地 降雨 集中 在 夏秋 季 , 降 雨量 的 增 
加 使 得 脲酶 活性 增加 ” ,这 使 得 土壤 细菌 网 络 中 与 
脲酶 联系 密切 的 部 分 细菌 更 加 活跃 ,土壤 细菌 在 脲 
酶 催化 作用 下 能 够 更 高 效 的 分 解 枯 落 物 ” ,而 枯 落 
物 的 分 解 能 够 显著 提高 土壤 过 氧化 氧 酶 活性 ,这 又 


使 得 土壤 细菌 网 络 中 与 过 氧化 氨 酶 联系 密切 的 部 分 
细 戎 更 加 活跃 ,过 氧化 氢 酶 能 酶 促 过 氧化 氨 分 解 为 
水 和 氧 , 从 而 解除 过 氧化 氨 对 植物 的 毒害 作用 
另外 , 桔 落 物 分解 过 程 中 转化 酶 在 葡萄 糖 转化 上 有 
着 突出 贡献 ” 。 土 壤 酶 活性 的 提高 有 利于 使 更 多 
土壤 细菌 通过 网 络 产 生 联系 ,共同 参与 樟 子 松林 物 


(1) M25 a 到 43 a, 呼 伦 贝 尔 沙 地 樟 子 松 人 工 
林 土 壤 细 菌 网 络 愈加 复杂 和 紧密 ,但 土壤 细菌 群落 
对 环境 扰动 的 敏感 性 增加 ,43 a 人 工 林 最 易 受 到 影 
响 。 土 壤 深 度 从 0~10 cm 到 10~20 em, 人工 林 土壤 
细菌 网 络 复杂 性 和 紧密 度 降低 ,土壤 细菌 群落 对 环 
境 扰动 的 敏感 性 下 降 ;草地 土壤 细菌 网 络 愈加 复杂 
和 紧密 ,土壤 细菌 群落 对 环境 扰动 的 敏感 性 增加 。 
与 草地 相 比 ,人 工 林 土壤 细菌 网 络 复杂 性 较 低 。 

(2) 呼伦贝尔 沙 地 不 同 林 龄 樟 子 松 入 工 林 土壤 
细菌 网 络 关键 菌 种 类 型 存在 差异 。25 a .34a 和 43 a 
分 别 识别 出 3 种 .1 种 .6 种 不 同类 型 的 关键 菌 种 ,其 
中 隶属 于 酸 杆 菌 门 的 关键 菌 种 数量 最 多 。 草 地 关 
菌 种 识别 出 1 种 ,隶属 于 匿 杆 菌 门 。 

(3) 呼伦贝尔 沙 地 樟 子 松 人 工 林 土 壤 细菌 网 络 
受 土壤 含水 量 SO ALI ER AA MEW 
Tc D aE eB A ak SG AU SE E RAE Eo 
酶 活性 驱动 , 且 土 壤 有 机 质 对 土壤 细菌 网 络 中 具有 
高 中 介 中 心性 的 细菌 影响 最 大 。 
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Soil bacterial networks in Pinus sylvestris var. mongolica plantations 
of the Hulunbuir Desert 


ZHANG Songan, LIU Xuan, ZHAO Peishan, GAO Guanglei™, ZHANG Ying", 
DING Guodong"", LIU Ye' REN Yue 
(1. School of Soil and Water Conservation, Beijing Forestry University, Engineering Research Center of Forestry 
Ecological Engineering, Ministry of Education, Beijing 100083, China; 2. Yanchi Ecology Research Station of the 
Mu Us Desert, Yanchi 751500, Ningxia, China; 3. Key Laboratory of State Forestry and Grassland Administration on 
Soil and Water Conservation, Beijing 100083, China) 


Abstract: To elucidate soil bacterial network interactions within Pinus sylvestris var. mongolica plantations in 
the Hulunbuir Desert. P. sylvestris plantations representing three different age groups (25 a, 34 a, and 43 a) and 
two soil layers (0- 10 and 10-20 cm) were selected to assess their soil bacterial networks using molecular 
ecological network analysis and data from a referenced grassland. The numbers of network edges increased, the 
average path length reduced from 25 a to 43 a. While the number of network edges reduced and the average path 
length increased with soil layers from 0-10 cm to 10-20 cm. Compared with the grassland, the P. sylvestris 
plantations had a lower network edge, and the soil bacterial network was less complicated. The soil bacteria were 
found to belong to the Acidimicrobiales, RB4/, and MB-A2- 108 in the 25 a plantation, Gaiellales in the 34 a 
plantation, and Gaiellales, RB4/, Subgroup 7, Subgroup 6, and DAI0I soil group in the 43 a plantation, 
Latescibacteria in the grassland. The soil bacterial network was significantly positively correlated with ammonia 
nitrogen, total nitrogen, microbial carbon content, and urease activities and significantly negatively correlated 
with invertase and catalase activities, soil water content and available phosphorus (P « 0.05). The soil organic 
matter had both positive and negative effects (P « 0.05). The soil bacterial network complexity and compactness 
increased from 25 a to 43 a. The opposite was found for the soil layers from 0-10 cm to 10-20 cm. Compared 
with the grassland, the soil bacterial network of the P. sylvestris plantation was less complicated. The keystone 
soil bacteria species were different among the three stand ages, and were greatest in stand 43 a. There were more 
keystone species in the plantation than the grassland. Soil physicochemical properties and enzymatic activity 
derived the soil bacterial network, and soil organic matter was the major influencing factor. This improved 
information contributed to a deep understanding of the soil bacterial community and provided a scientific and 
technological basis for the sustainable management of P. sylvestris plantations in the Hulunbuir Desert. 
Keywords: soil microbes; molecular ecological network; stand age; soil depth; Pinus sylvestris var. mongolica; 
the Hulunbuir Desert 


